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Abstract 



The QCD vacuum condensates in the Operator Product Expansion are extracted 
from the final ALEPH data on vector and axial-vector spectral functions from 
r-decay. Weighted Finite Energy Sum Rules are employed in the framework 
of both Fixed Order and Contour Improved Perturbation Theory. An overall 
consistent picture satisfying chirality constraints can be achieved only for values 
of the QCD scale below some critical value A ~ 350 MeV. For larger values of 
A, perturbation theory overwhelms the power corrections. A strong correlation 
is then found between A and the resulting values of the condensates. Reason- 
able accuracy is obtained up to dimension d = 8, beyond which no meaningful 
extraction is possible. 
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1 Introduction 



In the absence of exact analytical solutions to QCD at the Fermi scale, approxi- 
mation methods are required to confront the theory with hadronic information. 
A particularly appealing method is that of QCD sum rules [I], based on the 
Operator Product Expansion (OPE) of current correlators, and on the notion 
of quark-hadron duality. Dispersion relations relate hadronic parameters, such 
as masses and couplings, to universal QCD parameters, e.g. the strong coupling 
constant in the perturbative QCD sector, and vacuum expectation values of 
gauge invariant quark and gluon field operators in the non-perturbative sector. 
These vacuum condensates are organized according to dimension, and their nu- 
merical values can be determined from data on hadronic spectral functions, e.g. 
as obtained from e^e~ annihilation, or r- lepton decay data. The former de- 
termines the vacuum condensates in the vector channel, while the latter allows 
for the extraction of both vector and axial-vector channel condensates. Given 
the important role they play in applications of QCD sum rules, it is essential 
to determine the condensates with reasonable accuracy. Experimental data on 
hadronic spectral functions reaches up to some maximum energy, beyond which 
one expects perturbative QCD (PQCD) to provide a reasonable description 
of the current correlator. This so called continuum threshold, sq, is typically 
So — 1 — 3 GeV^. The vacuum condensates are then expressed in terms of a 
dispersive integral of the experimental spectral function up to s = sq, followed 
by a calculable integral of the PQCD expression of the correlator. The latter 
obviously involves the strong coupling constant as{s). Hence, there results a 
correlation between the extracted values of the vacuum condensates and the 
input value of the QCD scale A entering as{s). The resulting sum rules based 
on Cauchy's Theorem are called Finite Energy Sum Rules (FESR). A valuable 
feature of this approach to the determination of the condensates is that, to first 
order in as{s), radiative corrections to the condensates do not induce mixing 
of condensates of different dimension in a given FESR [2j. This decoupling is 
absent in other sum rules, e.g. Laplace or Gaussian, which then complicates 
substantially the analysis. 

Some of the early determinations of the condensates [3] used values of A in 
the range A ~ 100 — 200 MeV, which is substantially below the current range 
A ~ 300 — 400 MeV [1]. Such high values of A seriously complicate the task 
of extracting the condensates, as we pointed out a few times in the past [S]-[3]. 
The problem is that with increasing A, PQCD becomes the dominant contri- 
bution to the FESR, and overwhelms the power corrections above some critical 
value Ac. We have estimated |5]-[6] this to happen for Ac ~ 330 — 350 MeV. 
Another problem affecting the extraction of the vacuum condensates is the slow 
saturation of dispersive integrals by the hadronic spectral function data. Un- 
like the situation in deep inelastic electron-proton scattering, where precocious 
scaling is observed, the approach to asymptotia in tau-lepton decay does not 
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appear to be as fast. The final ALEPH [7] data for the chiral spectral function 
v{s) — a{s) shows clearly that this function has not yet reached its asymptotic 
form dictated by PQCD, i.e. it does not vanish, even at the highest energies 
attainable in r-decay. If the asymptotic regime had been reached precociously, 
let us say at s ~ 2 GeV^, then it would have been straightforward to calculate 
the non-perturbative condensates with the help of the Cauchy integral. Since 
this is not the case, some method to improve convergence must be used. We 
have shown [S], [H]-[S] that in the framework of FESR this can be done after 
suitably reducing the impact of the high energy region in the dispersive integral 
by using weighted sum rules. We used the data in a weighted linear combination 
of the first two Weinberg sum rules, which follows from the absence of dimension 
d = 2 and d — A operators in the chiral correlator, to demonstrate the preco- 
cious saturation of the FESR and the remarkable effectiveness of the method. 
Motivated by this success, we determined a number of QCD condensates by 
making maximal use of the absence of d = 2 and d = 4 chiral operators, and 
requiring strong stability, i.e. we varied the radius sg in the Cauchy integral 
beginning at the end of r-decay phase space and required that the condensates 
calculated from the data should be reasonably constant for all sq, in some finite 
region, including the end of phase space [3] . The results are in agreement with 
most independent determinations |10| . 

The next step is to proceed from chiral condensates to the extraction of in- 
dividual vector and axial-vector condensates. In this task, it is important to 
achieve a consistent picture, e.g. the d = 4 condensates should be chiral sym- 
metric (neglecting the small term niq < qq >) a.s they are proportional to the 
gluon condensate. The same symmetry is expected of potential d = 2 terms 
in the OPE, if present at all [S], [TT]. In addition, individual values for vector 
and axial-vector condensates of dimension d — 6 should be consistent with in- 
dependent determinations of chiral condensates of the same dimension [9]-|10|. 
This already poses a problem, as the chiral sum rules are independent of A (for 
vanishing light quark masses), while the individual FESR are expected to intro- 
duce an appreciable correlation between A and the condensates. The question 
is how strong is this correlation. In this context, an issue often overlooked is 
the manifest difference between the saturation of FESR by data in the vector 
as opposed to the axial-vector channel. The former involves the rho-meson res- 
onance peak at relatively low energy with a narrow width, while the latter has 
the very broad oi resonance above 1 GeV. As a result of this, the difference 
between the hadronic integral and the PQCD contribution, which essentially 
gives the condensates, is saturated differently in the two channels. For instance, 
for dimension d — 2, saturation in the vector channel is from above, while in 
the axial-vector case it is from below. An important additional constraint is 
provided by data on the r hadronic width R^-, which involves QCD perturba- 
tion theory as well as vacuum condensates of dimension d = 6 and d = 8. The 
theoretical expression for Rt involves a different integration kernel as the FESR 
giving the condensates; it can then be considered as a relatively independent 
constraint. Developments in the perturbative QCD sector, i.e. Contour Im- 
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proved Perturbation Theory (CIPT)[T2] should also be incorporated into the 
analysis. In fact, GIFT has been shown to provide a faster approach to asymp- 
totia than Fixed Order Ferturbation Theory (FOFT) . 

In this paper we carry out this self consistent analysis using the latest and final 
ALEFH data [7 on the separate vector and axial vector spectral functions. This 
data differs from earlier ALEFH versions in that all data (as opposed to partial 
data sets) collected from 1991-1995 was used, and in that an increase in statistics 
has resulted in a reduction of the overall errors. It is important to mention at 
this point that using previous ALEFH data sets we were unable to obtain a con- 
sistent picture of the dimension d = A condensates in the framework of FOFT. 
In fact, their signs differed, while they should have been identical. This problem 
is not present if the latest ALEFH data is used. In Section 2 we use FOFT, 
and in Section 3 we determine the condensates in the framework of GIFT. We 
favour a separate determination for each dimension and for each channel, for 
a particular value of A, as opposed to a single massive x-squared fit based on 
some exhaustive set of sum rules. The latter poses the danger of masking the 
correlation between perturbative and non-perturbative contributions to the dis- 
persive integrals. At the same time, it is very likely unstable with respect to the 
introduction of additional higher dimensional condensates. It can also overlook 
chirality constraints, necessary to achieve an overall consistent picture. As in 
the chiral case we find that a consistent picture can only be achieved on the 
basis of weighted, often called pinched, sum rules. An essential input to such 
sum rules is the absence of dimension d — 2 condensates. We shall show that 
this is indeed a viable assumption. The conclusions are presented in Section 4. 



2 Fixed Order Perturbation Theory 

We begin by defining the vector and axial-vector correlators 




: J d^xe^"^ <0|T(V;,(a;) 




(1) 




: j d^xe'-i^ <0|T(A^(x) 4(0))|0> 



(2) 
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where =: {wyf^u- d'y^d) : /2, and =: (wyf^j^u- d'j^jsd) : /2. Considering 
these (charge neutral) currents implies, in perturbative QCD, the normalization 

i ImUr^ is) = i Imn^" - ^ (l + - + ■••) ■ (3) 
The OPE of the current correlators Eqa.©-© may be written as 



N=0 ^ 

where the term N — Q stands for the purely perturbative QCD expression. An 
alternative definition of the condensates, used often in the literature is 

C2jv(g',M') < 02Jv(/i') > =87r2 d2Ar(Q')|v,A. (5) 
In FOPT one first calculates the Cauchy contour integral for fixed as{^) 

ds Tl{s)\^^/^ = f '° ds - Im nis)\vA , (6) 
Jo I" 

after which one performs the Renormalization Group (RG) improvement. This 
leads to 



fsn 1 N+1 

{- fC2N+2 < 62JV+2 87r2 / ds s^-Im n(s) - . ° , , /^(sq) , (7) 

Jo + 1) 

where iV = 0, 1, 2, the label V,A is omitted in the sequel, and In{so) is given 
by the left hand side of Eq.®. An explicit calculation to five- loop order in 
PQCD yields 
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Figure 1: The channel-averaged dimension d = 2 vacuum condensate from 
Eq.Q for (a): A = 300 MeV, and (b): A = 350 MeV. 
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where for three flavours, as used throughout in this paper, and in the AIS 
scheme, the coefHcients are [13]: /3i = -|, /^a = /^s = -^3 = 1-6398, 

-F4 = —10.2839. The strong couphng constant in Eq.®, as(so) = cts{so)/'!T, to 
five-foop order is [Hj 



(4)/ ^ 
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with L = hi(so/A^), and 
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with Ca — 1.202. The five- loop constant is not yet known. We have estimated 
it [5] assuming a geometric series behaviour for those constants not determined 
by the renormahzation group, i.e. ~ /e|/A:i ~ 25, with ki = F3 and /c2 = 
F4 + tt'^ 01/12. This is in good quahtative agreement with other estimates |15| . 

In order to compute the hadronic integral in Eq.(7) we have used the latest 
and final ALEPH data set [7]- This set includes together with each data point 
the statistical error, and also the full error correlation matrix. The integration 
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interval is divided in bins, and the hadronic integral is computed as a function 
of So in the standard fashion by addition, taking into account error propagation 
and correlation. From Eq.([7|) with = 0, we determined the dimension d = 2 
condensate in the vector and the axial-vector channel, for A = 300 — 350 MeV. 
The results do not yield a consistent picture, as the vector channel condensate 
is small and negative in the stability region around sq — 1.7 — 2.3 GeV^, while 
in the axial-vector channel it is mostly positive and numerically much larger. 
However, constraining this condensate to be chiral-symmetric, and taking the 
average in the two channels leads to a very stable result, consistent with zero, 
as shown in Fig. 1. The criterion of stability, which is standard in this kind of 
analysis, is quite simple. Quantities such as the vacuum condensates should be 
indepenedent of sq. Therefore, their values are read in the region where they 
show the least dependence on sq, which should be high enough to allow for the 
onset of perturbative QCD. Numerically the result is 



^ " " . — A , ^ 0.05 ± 0.05 GeV^ . (12) 
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C2 < O2 > k + C2 < O2 > 



For higher dimensional condensates, d > 4, results from these standard FESR 
are inconclusive. For instance, the d = 4 gluon condensate is not positive definite 
in either channel. This is an indication that the asymptotic QCD regime has 
not been reached at the scales accessible in r-decay, at least in the vicinity 
of the positive real axis of the contour integral. An attempt has been made 
to improve this situation by introducing an integration kernel in Eq.([5]) that 
vanishes at s = sq, i.e. one makes the substitution ^ [1 — (s/sq)^] for 
A'' = 1, 2, .... A meaningful use of these so called pinched moments relies heavily 
on the assumption that there is no dimension d — 2 condensate. Results from 
these weighted FESR are much better, provided A ~ 300 MeV. For larger values 
of A the condensates are consistent with zero. Numerically we obtain 



C4 < O4 > |y,A = 0.10 ±0.05 GeV^ (A = 300 MeV) , (13) 
Ce<Oe>\v = -0.5 ± 0.2 GeV^ (A = 300 MeV) , (14) 

Ce < Oe > U = -0.3 ± 0.2 GeV^ (A = 300 MeV) . (15) 

In the case of d = 8, only the vector channel condensate shows reasonable 
stability and has the value 
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Figure 2: The dimension d = A vacuum condensate from the vector channel in 
CIPT, Eq.(l22l), for (a): A = 300 MeV, and (b): A = 350 MeV. 



0.6 ± 0.2 GeV* 



(A= 300 MeV) 



(16) 



The graphical dependence of the condensates on the continuum threshold sqj 
being qualitatively similar in FOPT as in CIPT, will be shown in the next 
section. There, we discuss the use of CIPT together with pinched moments 
which does lead to a more consistent determination of the first few vacuum 
condensates with somewhat improved accuracy. 



3 Contour Improved Perturbation Theory 



Unlike the case of FOPT, where as(so) is frozen in Cauchy's contour integral, 
and the RG is implemented after integration, in CIPT ag is running and the 
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Figure 3: The dimension d = A vacuum condensate from the axial- vector channel 
in CIPT, Eq.jSSI), for (a): A = 300 MeV, and(b):A = 350 MeV. 

RG is used before integrating. It has been shown that in r-decay CIPT is 
superior to FOPT. To proceed we introduce the Adler function 

D{s) ^ -s-^n{s) . (17) 

From Cauchy's theorem, and after an integration by parts, a relation between 
the contour integrals of n(s) and D{s) is easily obtained 



/ ds ( A)A^ n(.) = ^ ^ / ^ _ ,s„^+i) Dis) . (18) 

After RG improvement, the perturbative expansion of the Adlcr function can 
be written as 
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Figure 4: The dimension d = 6 vacuum condensate from the vector channel in 
CIPT, Eq.(l22l), for (a): A = 300 MeV, and (b): A = 350 MeV. 



as(-s) 



(19) 



m=0 



where [16] iiTo = -^^i = 1, ^2 = -F3 = 1-6398 ,^^3 = ^4 + ^Pf = 6.3710,, for 
three flavours, and the unknown coefficient = ^03 ~ 25. Next, we introduce 
the moments 



Mat (so) 
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Km [lN + l,niiso) - Io,rn{so)] , (20) 



where 
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as{-s). 



LN.m =1 f ds ( — ) 

J\s\=so So ' TT 

The vacuum condensates can then be expressed as 



(21) 
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Figure 5: The dimension d = 6 vacuum condensate from the axial- vector channel 
in CIPT, Eq.(l22l), for (a): A = 300 MeV, and (b): A = 350 MeV. 
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where a pinched integration kernel has been introduced. The l.h.s. of Eg. ((22)) 
includes radiative corrections up to order 0{as). Mixing of operators of dimen- 
sion lower or higher than the leading term occurs only at order 0{a1) [2], which 
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can be safely neglected. In the framework of CIPT [T^ the contour integrals, 
Eq. (j2ip . are calculated numerically using for as a numerical solution to the 
RG equation 



f^^^^^f3n[as{f^r'. (23) 

n—l 



We have used a Modified Euler method to solve the RG equation at each inte- 
gration point in the complex s-plane, and a single step numerical integration to 
compute the contour integral. 

Using Eq. (j22p we first determine the dimension d — 2 vacuum condensates in 
the vector and in the axial- vector channel. Results are comparable to those from 
FOPT, so the numerical value of the channel-average condensate is still that of 
Eq. p^ . Next, for d — 4, and assuming C2 < O2 > \v,A — 0, the results in the 
vector and the axial-vector channel are shown in Figs. 2 and 3, respectively, for 
the two values A = 300 MeV, curve (a), and A = 350 MeV, curve (b). Although 
results in the axial-vector channel become consistent with those in the vector 
channel only close to the end-point, at least such a consistency is achieved. 
Numerically we find 



„ „ I _j 0.15 ± 0.04 GeV* (A = 300 McV) 
Ci<Oi>\v.A-i 0.07 ± 0.02 GcV* (A = 350 McV) ' ^^^^ 



For larger values of A the process breaks down as C4 < O4 > |y vanishes and 
C4 < O4 > \a becomes negative. This confirms our earlier claims about the 
existence of a critical value Ac ~ 330 — 350 MeV beyond which no meaningful 
determination of the condensates is possible [S]-[3]. Proceeding to dimension 
d — 6, results are shown in Figs. 4 and 5, for the same two choices of A. The 
condensate in the vector channel is reasonably stable giving the values 



„ „ I -0.63 ± 0.12 GeV^ (A = 300 MeV) 

C6<06>k-i _o.35±0.12GeV« (A = 350 MeV) ' 



Results in the axial- vector channel change sign at sq ~ 2.1 — 2.2 GeV^, and are 
then barely conclusive. However, we can use the reasonably accurate determi- 
nation [9] of the chiral condensate of d = 6 



Ce < Oe > k-A = -0.18 ±0.04, (26) 
together with Eq. ([25| to obtain 
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Figure 6: The dimension d = 8 vacuum condensate from the vector channel in 
CIPT, Eq.(l22l), for (a): A = 300 MeV, and (b): A = 350 MeV. 



-0.45 ±0.13 GeV^ 
-0.17 ±0.12 GeV^ 



(A = 300 MeV) 
(A = 350 MeV) 



(27) 



These results are in remarkable agreement with the values of Cg < Oq > \a, if 
one were to read off the condensate close to the end point (see Fig. 5) 



Ce <Oe>\A = 



-0.50 ±0.15 GeV'' 
-0.18 ±0.13 GeV^ 



(A = 300 MeV) 
(A = 350 MeV) 



(28) 



However, it should be kept in mind that the result Eq. (|26p is independent of 
A, as there is no perturbative contribution to chiral sum rules (for vanishing 
light quark masses). The sign of Ce < Oe > U as obtained above does not 
agree with the so called Vacuum Saturation (V.S.) approximation to estimate 
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the four-quark contribution to Cg < Oe >• In fact, according to V.S. the ratio 
of vector to axial- vector condensates is negative 



Ce < Oe > \v ^ T_ .^q^ 

Cr <Ofi > u 11 ■ ^ ' 



Reading off Cg < Oe > U at lower values of Sqi where it is positive, would give 
Ce < Oe > U =i 0.2 - 0.3 GeV^ leading to a chiral value -(0.65 - 0.85) GeV*^ 
in serious disagreement with Ea. (j26p . From Fig. 5 one could always speculate 
that the fully asymptotic region may have not been reached at 3 GeV^, and 
that a sign change in that region could take place. It should be mentioned that 
using the previous versions of the ALEPH data base, we obtain a positive sign 
for Oe < Oe > |a, and a negative sign for Oe < Oe > |y, using FOPT and a 
pinched moment. However, this data base leads to different signs for C4 < O4 > 
in the vector and axial-vector channels. 



Not surprisingly, the situation with the dimension d = 8 condensate is consid- 
erably worse. Results in the axial-vector channel are totally inconclusive: the 
condensate changes sign in the middle of the integration region for both values 
of A. For instance, for sq > 1.8 GeV^, Os < Og > \a rises monotonically from 
zero to about 2.5 GeV^ for A = 300 MeV. The vector channel provides an 
acceptable result only for A = 300 MeV, as shown in Fig. 6. Reading values in 
the stability region so — 1.5 — 2.4 GeV^ gives 



Og < Og > |y = 0.8 ± 0.2 GeV® (A = 300 MeV) . (30) 

In this case it is not possible to use this result, together with the independent 
determination of the condensate Og < Og > \v-a in the chiral limit, in order to 
obtain Og < Og > \a- In fact, at d = 8 there are chiral symmetric contributions, 
proportional to the gluon condensate and to the single quark condensate, which 
vanish in the chiral condensate J.7j . 

We finally turn to the r-ratio in either the vector or the axial-vector channel 



EW 



ds 
So 



l-3(-)2 + 2(-)3 
So So 



-Imn{s) . (31) 



The theoretical expression of which follows from Eqs.(|l]) and (|^D|) . again for 
either the vector or the axial-vector channel, is 
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Figure 7: The ratio in the vector channel. Data points are from Eq. ((3T|) . 
Theoretical curves are from CIPT, Eq.^, for A 300 MeV, and (a): 
Ce < Oe >= Cs < Os >= 0, (b): Ce < Oe > \v = -0-45 GeV^ and 
Cg < Os > |y = 0.78 GeV^ 



The results for R.^ in the vector channel are shown in FiglT] Curve (a) cor- 
responds to Eq.(I32]) for Cg < Oe >= Cs < Os >= 0, and curve (b) for 
Ce < Oe >= -0.45 GeV^ Cg < Og >= 0.78 GeV^ The latter values are in 
very good agreement with the results Eas. (|25)) and ([30]). This is an independent 
confirmation of the strong correlation between the values of the condensates and 
the value of A. The situation in the axial-vector channel is inconclusive partly 





(32) 
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because of the large uncertainties in the condensates, and partly because the 
data approach the asymptotic regime from below, rather from above as in the 
vector channel. 



4 Conclusion 



In the framework of FOPT, and using standard FESR without the pinched 
kernel [1 — (s/sq)^], the dimension d — 2 condensate is consistent with zero 
after requiring it to be chiral symmetric, and averaging over the two channels. 
However, the result, Eq. (fT2)) . does not exclude a very small value for this conden- 
sate. For dimension d > 4 no self consistent picture is achieved. For instance, 
for d = 4 the gluon condensate comes out with different signs in the vector 
and axial- vector channels. This overall situation improves with the inclusion 
of the pinched kernel in the FESR, but then only if A ~ 300 MeV. Larger 
values of A yield condensates consistent with zero. This agrees with our earlier 
observations about the existence of a critical value Ac ~ 330 — 350 MeV beyond 
which no meaningful determination of the condensates is possible [5]-[6j. In the 
framework of CIPT, and using weighted FESR, the overall picture improves. 
It is possible to accommodate larger values of A up to A ~ 350 MeV, and to 
determine the condensates in the vector and axial-vector channels with reason- 
able accuracy. It is also possible to achieve consistency in terms of chirality: 
the gluon condensate turns out to have the same sign in both channels, and the 
difference between the d = 6 condensates is consistent with independent deter- 
minations from chiral sum rules. The sign of Cg < Oe > U, though, disagrees 
with the Vacuum Saturation approximation, provided Cg < Oe > |v is negative, 
as we have found here. Of course, a change of sign in Cg < Oe > U beyond 
3 GeV^ cannot be excluded. For dimension d = 8, only the vector condensate 
can be determined unambiguously, and then only for A = 300 MeV. Given the 
strong correlation between A and the condensates, we wish to suggest that in 
applications of QCD sum rules, values for the latter be used in conjunction with 
the corresponding values of A. Because of this correlation, we find little merit 
in averaging the results obtained from this analysis at different values of A. 



References 

[1] For a recent review see e.g. P. Colangelo, A. Khodjamirian, in: "At the 
Frontier of Particle Physics/ Handbook of QCD'", M. Shifman, ed. (World 
Scientific, Singapore 2001), Vol. 3, 1495-1576. 

[2] G. Launcr, Z. Physik, C 32 (1986) 557. 



16 



[3] R.A. Bertlmann et al, Z. Phys. C 39 (1988) 231; C.A. Dominguez and J. 
Sola, Z. Phys. C 40 (1988) 63; V.Gimcncz, J.A. Pefiarrocha, J. Bordes, 
Phys. Lett. B 223 (1989) 245. 

[4] Particle Data Group, W. -M. Yao, J. Phys. G 33 (2006) 1. 

[5] C.A. Dominguez, K. Schilcher, Phys.Rev. D 61 (2000) 114020. 

[6] C.A. Dominguez, K. Schilcher, Phys. Lett.B 581 (2004) 193. 

[7] ALEPH Collaboration, S.Schael et ai, Phys. Kept. 421 (2005) 191. 

[8] C.A. Dominguez, K. Schilcher, Phys. Lett.B 448 (1999) 93. 

[9] J. Bordes, C.A. Dominguez, J. Peiiarrocha, K. Schilcher, J. High Energy 
Phys. 02 (2006) 037. 

[10] B.L. loffe, Progr. Part. Nucl. Phys. 56 (2006) 232; S. Narison, Phys. Lett.B 
624 (2005) 223; J. Rojo, J.I. Latorre, J.High Energy Phys. 0401 (2004) 
055; V. CirigHano, E. Golowich and K. Maltman, Phys.Rev. D 68 (2003) 
054013; B.V. Geshkenbein, Phys. Rev. D 67 (2003) 074006; M. Davier, A. 
Hocker and Z. Zhang, hep-ph/0507078 (2005). 

[11] S. Narison, Phys. Lett.B 300 (1993) 293; C.A. Dominguez, Phys. Lett.B 
345 (1995) 291. 

[12] A.A. Pivovarov, Sov. J. Nucl. Phys. 54 (1991) 676; F. Le Diberder and A. 
Rich, Phys. Lett. B 286 (1992) 147; M. Jamin, J. High Energy Phys. 09 
(2005) 058. 

[13] T. van Ritbergen, J.A.M. Vermaseren, S.A. Larin, Phys. Lett. B 400 (1997) 
379. 

[14] K.G. Chetyrkin, B.A. Kniehl, M.Steinhauser, Phys. Rev. Lett. 79 (1997) 
2184. 

[15] F. Le Diberder and A. Rich, Phys. Lett. B 286 (1992) 147; J. Ellis, M. 
Karhner and M.A. Samuel, Phys. Lett. B 400 (1997) 176; A.L. Kataev and 
V.V. Starshenko, Mod. Phys. Lett. A 10 (1995) 235; P.A. Baikov, K.G. 
Chetyrkin and J.H. Kuhn, Phys. Rev. D 67 (2003) 074026. 

[16] K.G. Chetyrkin, A.L. Kataev and F.V. Tkachov, Phys. Lett. B 85 (1979) 
277; M. Dine and J. Sapirstein, Phys. Rev. Lett. 43 (1979) 668; W. Celmas- 
ter and R. Gonsalves, Phys. Rev. Lett. 44 (1980) 560; S.G. Gorishny, A.L. 
Kataev, and S.A. Larin, Phys. Lett. B 259 (1991) 144; L.R. Surguladze 
and M. Samuel, Phys. Rev. Lett. 66 (1991) 560. 

[17] D.J Broadhurst and S.C. Generalis, Phys. Lett. B 165 (1985) 175. 



17 



